Search for: All records

We forward-model mass-weighted stellar ages (MWAs) and quiescent fractions (fQ) in projected phase space (PPS), using data from the Sloan Digital Sky Survey, to jointly constrain an infall quenching model for galaxies in log (Mvir/M⊙) > 14 galaxy clusters at z ∼ 0. We find the average deviation in MWA from the MWA–M⋆ relation depends on position in PPS, with a maximum difference between the inner cluster and infalling interloper galaxies of ∼1 Gyr. Our model employs infall information from N-body simulations and stochastic star-formation histories from the universemachine model. We find total quenching times of tQ = 3.7 ± 0.4 Gyr and tQ = 4.0 ± 0.2 Gyr after first pericentre, for 9 < log (M⋆/M⊙) < 10 and 10 < log (M⋆/M⊙) < 10.5 galaxies, respectively. By using MWAs, we break the degeneracy in time of quenching onset and time-scale of star formation rate (SFR) decline. We find that time of quenching onset relative to pericentre is $t_{\mathrm{delay}}=3.5^{+0.6}_{-0.9}$ Gyr and $t_{\mathrm{delay}}=-0.3^{+0.8}_{-1.0}$ Gyr for 9 < log (M⋆/M⊙) < 10 and 10 < log (M⋆/M⊙) < 10.5 galaxies, respectively, and exponential SFR suppression time-scales are τenv ≤ 1.0 Gyr for 9 < log (M⋆/M⊙) < 10 galaxies and τenv ∼ 2.3 Gyr for 10 < log (M⋆/M⊙) < 10.5 galaxies. Stochastic star formation histories remove the need for rapid infall quenching to maintain the bimodality in the SFR of cluster galaxies; the depth of the green valley prefers quenching onsets close to first pericentre and a longer quenching envelope, in slight tension with the MWA-driven results. Taken together these results suggest that quenching begins close to, or just after pericentre, but the time-scale for quenching to be fully complete is much longer and therefore ram-pressure stripping is not complete on first pericentric passage.

We explore models of massive (>1010 M⊙) satellite quenching in massive clusters at z ≳ 1 using an MCMC framework, focusing on two primary parameters: Rquench (the host-centric radius at which quenching begins) and τquench (the time-scale upon which a satellite quenches after crossing Rquench). Our MCMC analysis shows two local maxima in the 1D posterior probability distribution of Rquench at approximately 0.25 and 1.0 R200. Analysing four distinct solutions in the τquench–Rquench parameter space, nearly all of which yield quiescent fractions consistent with observational data from the GOGREEN survey, we investigate whether these solutions represent distinct quenching pathways and find that they can be separated between ‘starvation’ and ‘core quenching’ scenarios. The starvation pathway is characterized by quenching time-scales that are roughly consistent with the total cold gas (H2 + H i) depletion time-scale at intermediate z, while core quenching is characterized by satellites with relatively high line-of-sight velocities that quench on short time-scales (∼0.25 Gyr) after reaching the inner region of the cluster (<0.30 R200). Lastly, we break the degeneracy between these solutions by comparing the observed properties of transition galaxies from the GOGREEN survey. We conclude that only the ‘starvation’ pathway is consistent with the projected phase-space distribution and relative abundance of transition galaxies at z ∼ 1. However, we acknowledge that ram pressure might contribute as a secondary quenching mechanism.

We model satellite quenching at z ∼ 1 by combining 14 massive (1013.8 < Mhalo/M⊙ < 1015) clusters at 0.8 < z < 1.3 from the GOGREEN and GCLASS surveys with accretion histories of 56 redshift-matched analogues from the IllustrisTNG simulation. Our fiducial model, which is parametrized by the satellite quenching time-scale (τquench), accounts for quenching in our simulated satellite population both at the time of infall by using the observed coeval field quenched fraction and after infall by tuning τquench to reproduce the observed satellite quenched fraction versus stellar mass trend. This model successfully reproduces the observed satellite quenched fraction as a function of stellar mass (by construction), projected cluster-centric radius, and redshift and is consistent with the observed field and cluster stellar mass functions at z ∼ 1. We find that the satellite quenching time-scale is mass dependent, in conflict with some previous studies at low and intermediate redshift. Over the stellar mass range probed (M⋆ > 1010 M⊙), we find that the satellite quenching time-scale decreases with increasing satellite stellar mass from ∼1.6 Gyr at 1010 M⊙ to ∼0.6−1 Gyr at 1011 M⊙ and is roughly consistent with the total cold gas (HI + H2) depletion time-scales at intermediate z, suggesting that starvation may be the dominant driver of environmental quenching at z < 2. Finally, while environmental mechanisms are relatively efficient at quenching massive satellites, we find that the majority ($\sim 65{\!-\!}80{{\ \rm per\ cent}}$) of ultra-massive satellites (M⋆ > 1011 M⊙) are quenched prior to infall.

Recent observations have shown that the environmental quenching of galaxies at z ∼ 1 is qualitatively different to that in the local Universe. However, the physical origin of these differences has not yet been elucidated. In addition, while low-redshift comparisons between observed environmental trends and the predictions of cosmological hydrodynamical simulations are now routine, there have been relatively few comparisons at higher redshifts to date. Here we confront three state-of-the-art suites of simulations (BAHAMAS+MACSIS, EAGLE+Hydrangea, IllustrisTNG) with state-of-the-art observations of the field and cluster environments from the COSMOS/UltraVISTA and GOGREEN surveys, respectively, at z ∼ 1 to assess the realism of the simulations and gain insight into the evolution of environmental quenching. We show that while the simulations generally reproduce the stellar content and the stellar mass functions of quiescent and star-forming galaxies in the field, all the simulations struggle to capture the observed quenching of satellites in the cluster environment, in that they are overly efficient at quenching low-mass satellites. Furthermore, two of the suites do not sufficiently quench the highest mass galaxies in clusters, perhaps a result of insufficient feedback from AGN. The origin of the discrepancy at low stellar masses ($M_* \lesssim 10^{10}$ M⊙), which is present in all the simulations in spite of large differences in resolution, feedback implementations, and hydrodynamical solvers, is unclear. The next generation of simulations, which will push to significantly higher resolution and also include explicit modelling of the cold interstellar medium, may help us to shed light on the low-mass tension.